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We discuss the fate of chiral symmetry in an extremely strong magnetic field B. We investigate 
not only quark fluctuations but also neutral meson effects. The former would enhance the chiral- 
symmetry breaking at finite B according to the Magnetic Catalysis, while the latter would suppress 
the chiral condensate once B exceeds the scale of the hadron structure. Using a chiral model we 
demonstrate how neutral mesons are subject to the dimensional reduction and the low dimensionality 
favors the chiral-symmetric phase. We point out that this effect, the Magnetic Inhibition, can be 
a feasible explanation for recent lattice-QCD data indicating the decreasing behavior of the chiral- 
restoration temperature with increasing B. 



PACS numbers: 11.30.Rd, 21.65. Qr, 12.38.-t 

Theorists have been pursuing the answer of questions 
in various extremes. What happens if the temperature 
is extremely high? According to the theory of the strong 
interaction, namely, quantum chromodynamics (QCD), 
chiral symmetry should be restored and color degrees of 
freedom should be released at the temperature T of the 
order of A QCD - 200 MeV [1-3]. If T is raised fur- 
ther, the electroweak phase transition should take place 
at T ~ 100 GeV. Along the same spirit many theorists 
are trying to clarify what could happen at extremely 
high baryon density. In the QCD asymptotic regime 
where the perturbative calculation should work, theo- 
retical considerations predict the color superconducting 
phases [2, 4, 5]. In particular the ground state should 
form the color-flavor locking (CFL) if the quark chem- 
ical potential, yU q , is sufficiently larger than the strange 
quark mass, while there may appear many other pairing 
patterns such as the uSC phase, the dSC phase, the 2SC 
phase, etc in the intermediate density region. It is still 
a big theoretical challenge to identify the correct phase 
structure of QCD on the whole /i q -T plane. 

Recently the QCD phase structure in the presence of 
strong magnetic field B has been revisited extensively for 
several good reasons [3]: First, QCD matter under strong 
B is worth thinking as a realistic situation in non-central 
collisions in the relativistic heavy-ion experiment as con- 
ducted at the Relativistic Heavy-Ion Collider (RHIC) and 
the Large Hadron Collider (LHC). A simple calculation 
gives us an order estimate for the produced magnetic field 
as eB > Aq CD at the RHIC energy [6]. Second, the B- 
effect is very similar to the baryon chemical potential in 
a sense that gluons have no direct coupling and the QCD 
equation of state is affected only through quark polariza- 
tion processes. The Monte-Carlo simulation with B is, in 
spite of the similarity to /i q , possible without the notori- 
ous sign problem [7, 8], which is a great theoretical ad- 
vantage. Third, it is certainly intriguing to address such 
a simple and well-defined question; what is the ground 



state of QCD matter when a very strong magnetic field, 
B > Aqcd, is applied? 

It should be tough in general to answer to such a 
question since the confmcment/deconfinement phenom- 
ena belong to the gluon dynamics and B-effects are then 
indirect. (See Ref. [9] for recent attempts.) This diffi- 
culty is common also in the finite- /i q analysis, which hin- 
ders the QCD phase diagram research also [10]. In other 
words, one could never reach the correct QCD phase dia- 
gram on the )J.q-T plane until one can establish a machin- 
ery to reveal the B-effects that should be under better 
theoretical control [11]. In a related context the new en- 
tanglement with finite B and n q is also an interesting 
research subject [12]. 

In contrast to the confinement sector, the properties 
of chiral symmetry reside in the quark part, and thus 
they are directly sensitive to the presence of B. It would 
be, therefore, more tractable to focus on chiral symme- 
try in the strong B limit. From this point of view of 
the interplay between chiral symmetry and B, the Mag- 
netic Catalysis is one of the most significant phenom- 
ena [13-15]: In chiral quark models such as the Nambu- 
Jona-Lasinio (NJL) model without the confinement ef- 
fect, chiral symmetry is spontaneously broken for a suf- 
ficiently strong coupling constant, i.e. G > G c . With 
the magnetic field, on the other hand, a non-zero value 
of the chiral condensate would be an inevitable conse- 
quence from the Landau zero-mode contribution regard- 
less of the value of G. Fermions are always massive in 
the presence of B, hence, even though they are massless 
at the Lagrangian level (i.e. the chiral limit). In fact, the 
(3+1) dimensional NJL model with U(1)l x U(1)r chiral 
symmetry [14] [this is the unbroken part of chiral sym- 
metry with B that breaks isospin symmetry explicitly] is 
defined by the Lagrangian, 
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where p = ^(d^+ieAp) with A» = (0,-yB/2,xB/2,0) 
(in the symmetric gauge) and B > 0. Then it has been 
established that the constituent quark mass is expressed 
analogously to the gap obtained in the BCS theory; 

2 eB ( 4tt 2 \ 

which is non-zero for any G. One can understand Eq. (2) 
in such a way that B is a catalyst to induce a non- 
vanishing chiral condensate, that is, the Magnetic Catal- 
ysis. In the NJL model as an effective description of 
QCD at low energy, G has the energy scale comparable 
to Aqcd, i-c G ~ Aqcd- This means, together with 
Eq. (2), that the -B-induced value of to 2 becomes appre- 
ciable for eB ~ Aq CD . 

It is natural, as suggested by the Magnetic Catalysis, 
that B should enhance the chiral-symmetry breaking, so 
that the critical temperature, T c , for chiral restoration 
should increase with increasing B, which is indeed the 
case in all chiral-model calculations (see Ref. [16] for a 
recent review). The latest lattice-QCD data, however, 
supports dropping behavior of T c as a function of B and 
there is no clear physical explanation for this. The main 
goal of the present work is to propose a new mechanism 
for chiral restoration at strong B that could be a feasible 
explanation for the lattice-QCD data. 

Our idea is as follows. Because U(1)l x U(1)r is spon- 
taneously broken down to U(l)v, the Nambu-Goldstone 
(NG) boson (i.e. ir°) must exist as a composite of 
fermions. It does not matter whether ir° is a tight bound- 
state particle of QCD or not, but here let us just call 
this NG boson ir° in our QCD-based convention. If B 
is extremely strong, fermions that form ir° are affected 
by B and eventually their motions are restricted along 
the B-direction. Hence, it should be conceivable that 7r° 
also undergoes the dimensional reduction to the (1+1)- 
dimensional dynamics. Once this happens, the sponta- 
neous chiral-symmetry breaking is prohibited according 
to Marmin- Wagner's theorem [17]. Such a possibility, 
that we name the Magnetic Inhibition, was considered 
partially in Ref. [14], but only the approximated form 
of the 7T° propagator was discussed there. In this work 
we will fully evaluate the ir° propagator to address its 
dynamical change in strong B and formulate the above- 
mentioned idea. 

For later convenience let us look closely at the 
derivation of Eq. (2) using the quark propagator on 
top of the vector potential A^ . The quark prop- 
agator is expressed as S(p) = J2^=o ^n(p)/{p 2 — 
to 2 ) with m 2 = m 2 + 2eBn and S n (p) = (jj>« + 
m)[P + A n {p 2 ± ) + P_A„_ 1 (p 2 L )] + ^B n (p 2 ± ). Here we 
have introduced several notations: p 2 = p\ — p 2 , p\ = 

pl+p 2 y , A n {p\) ee 2e- 2z (-l) n Ll 0) (4 2 ), and B n {p\) = 

4e- 2z (-l) n L^ 1 2 1 (4z) with z = p 2 L /{2eB) and the pro- 
jection operators; P± = (1 ± ij 1 j 2 )/2. The generalized 




FIG. 1. Schematic expression of the gap equation. The solid 
curve represents the quark loop with the quark propagator 
S(p) and the four-Fermi coupling G and the dashed curve 
represents the pion loop with the pion propagator D n and 
the coupling r CTvr7r . 



Laguerrc polynomials are defined as usual by L„ (x) = 
(e x x- a /n\)(d n /dx n ){e- x x n+a ) [14, 18]. 

The most important ingredient to study the Magnetic 
Catalysis is the gap equation at the quark one-loop level, 

which is regularized in the proper-time method. 

This represents the tadpole diagram in terms of quarks 
(see the loop diagram with a solid line in Fig. 1). Ob- 
viously m = is a solution of the above equation. We 
have to make a comparison in energy to locate the ground 
state that should have the lowest energy. It is a straight- 
forward exercise to evaluate the effective potential by in- 
tegrating the gap equation with respect to to, which leads 
to V(m) = to 2 /(2G) + Vq(m), where the quark part is 

VJm) = I", / "T e ' m ' S coth(eBs) + (const.) . (4) 
8tt 2 J 1/A 2 s 2 

The onset for the spontaneous chiral-symmetry breaking 
in the chiral limit is immediately located from the sign 
of the potential curvature, i.e. the coefficient of the m 2 - 
term in the potential. In the B = case, we can expand 
as T^q (to) ~ — (A 2 /87r 2 )m 2 + (const.), from which we can 
conclude that chiral symmetry is spontaneously broken 
only for GA 2 > 47r 2 . [Note that we did not consider 
the color factor here.] In the limit of eB 3> A 2 , on the 
other hand, only the lowest Landau level (Landau zero- 
mode) can contribute to the gap equation, or, we can 
approximate the gap equation as coth(eBs) ~ 1 to find 

pR / p 1 ^ A 2 \ 

V q (m) ~ |+r ( A 2 - to 2 In + 0(m 3 ) . (5) 

871-^ \ TO / 

We see that the potential curvature has a logarithmic sin- 
gularity at to = and thus the curvature can be always 
negative for sufficiently small m, which means that the 
symmetric state with to = cannot be realized. This is 
how the Magnetic Catalysis works. The extremal point 
of m 2 /(2G) + Vq(m) with Eq. (5) gives a gap equation 
whose solution reads to 2 = e -7 A 2 e~ 47r /( eSG ). We note 
that the difference in the overall coefficient from Eq. (2) 
originates from whether eB ^> A 2 or not. As long as 
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eB <C A 2 , the dynamical quark mass is characterized by 
eB, but once eB exceeds the order of A 2 , the quark mass 
squared is no longer proportional to eB but suppressed 
by another (smaller) scale of A 2 . 

Now let us proceed to the calculations including the 
pion-loop effects under strong B. Because charged pions 
are as massive as eB, we can simply discard tt^ and focus 
only on 7r°, which justifies the usage of our simple model 
setting with only U(1)l x U(1)r except for the color fac- 
tor. As long as eB is small as compared to the pion size, 
we can treat 7T° as a point particle as in the chiral per- 
turbation theory [15, 19]. However, ir° is a composite 
particle, and it is conceivable that the dispersion relation 
of 7T° should be significantly modified by B. We can con- 
cretely investigate this by constructing it dynamically 
in the present model. In the conventional random phase 
approximation the pion propagator is 



iD~ 1 (p)=~+i 



d A k 



ti[ 75 S(k) 75 S(p + k)] . (6) 
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After some (tedious) calculations we can find the follow- 
ing expression: 
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where we introduced new notations to indicate 
some combinations of the Lagucrrc polynomials, 
Fl (z) = [L { r l \z)f + (n/i)[ife (*)] 2 and F&z) = 
(z/l)[L\^ l+1 \z)} 2 + (n/z)[i ; ( "" i ~ 1) (z)] 2 . Also, we de- 
fined 



n2(pjj,m 
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(2tt) 2 k I -ml (fc||+pn) 5 



■p\\) — u H 

(8) 

2eBl. The ex- 
pression looks complicated and it would be convenient to 
approximate it as 



with 



2eBn and m 2 = m 2 



iDZ 



z: 



2 2 
V±P± 



(9) 



that is motivated from an expansion valid for p\ < eB 
and p 2 < m 2 . The explicit computation gives Z~ x = 
(l/8TT 2 )[eB/m 2 + \n{e~~< A 2 /2eB) - ^(1 + m 2 /2eB)} and 
v\ = (Z 7r /8TT 2 )\n(e-' f A 2 /m 2 ) with tp(x) being the 
digamma function [14]. We note that the physical pion 
mass is given by m 2 , = Z^ih 2 using the bare pion mass 
rh^ in Eq. (7). 

To see how Eq. (9) works, we numerically evaluate the 
full propagator (7) to make a plot for the dispersion rela- 
tion in Fig. 2. As long as Eq. (9) is a sensible approxima- 
tion of Eq. (7), the zero of the pion propagator inverse 
should behave like p| = v\p\ + m 2 , which is clearly con- 
firmed in Fig. 2. Furthermore, this type of the dispersion 
form persists even for p 2 ^ > eB and/or p 2 > m 2 . 




FIG. 2. Zero of the pion propagator inverse (7) as a function 
of pjj, p\, and eB in the unit of A 2 . We chose m = 0.5A to 
avoid unphysical threshold effects. The slope of p 2 against 
p\ corresponds to the transverse velocity v\. Clearly the 
expanded form (9) is a good approximation even for p\ > eB 



and/or p<< > 



Concerning the properties of Z^ and v\, here, the es- 
sential point is that Z~ x ~ 0(eB) and thus Z v goes 
smaller with increasing eB. This makes the transverse 
velocity of 7r° behave as v\ ~ Z^ ~ l/(eB) that goes 
smaller accordingly. Such vanishing behavior of v\ is 
nothing but the concrete realization of the dimensional 
reduction from the (3+1)- to the (l+l)-dimcnsional dy- 
namics. 

We would make a remark about the nature of the di- 
mensional reduction for quarks and pions. One might 
have thought, at a first glance, that the dimensional re- 
duction with v\ — > in Eq. (9) seems to be a different 
situation from quarks under strong B. In an intuitive 
picture quarks are trapped by B and the transverse mo- 
tion is highly restricted, so that quarks can move only 
along B, which is how the dimensional reduction occurs 
for quarks. On the other hand, the neutral pion costs 
no energy to move and thus travels freely in the trans- 
verse directions when v 2 ^ = 0. Such intuitive descrip- 
tions may sound far different but the underlying physics 
is common. Actually the Landau level for quarks is just 
a quantum number, and even for the Landau zero-mode 
for example, the pj^-integration should be carried out, 
which picks up the Landau degeneracy factor, eB/(2ir). 
In the same manner the p^-integration for pions should 
count the density of states. 

We shall explicitly go into the computation of the pion 
loop as depicted by the dashed line in Fig. 1. For pjj > 
4m 2 the full propagator given by Eq. (7) would suffer 
from the threshold effect associated with the tt° — > qq 
decay which is unphysical due to the lack of confinement. 
We can evade this artifact by keeping the approximate 
form (9), the validity of which is checked in Fig. 2. We 
also cut off the transverse momentum integration in the 
range p\ < eB, which is reminiscent of the coefficient in 
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Eq. (2); the transverse degeneracy factor should be either 
eB or A 2 that is smaller than the other. The microscopic 
origin of this cutoff by eB is the p^-dcpcndence in v\. 
Not to rely on model-dependent details, we postulate the 
form (9) and introduce a sharp cutoff with an unknown 
parameter £ as 



P±<SeB 



d 4 p 
(2n) 4 



r(p)Ar(p) 



2m 



2m 



p\<(,eB 



l/AJ 



d 4 p i 
(^O 1 p\ - v\ V \ 



d 4 p 



/ ~2 , 2 2 i 2 \ 



(10) 



where the triple meson vertex is given by r CT7r7r (p) = 
— (6 1 'Sm)iD~ 1 (p) that we approximate at vanishing mo- 
mentum by r (T7r7r (0) = 4m/Z 7r . We have used the Wick 
rotation from p to Euclidean p and implemented the 
proper-time regularization again with the UV cutoff A w 
that is in principle related to the cutoff A in the quark 
sector. This p-integration is finite and results in the fol- 
lowing expression; 



id 



8n 2 v 2 ± 



( 1 _ e -»feBt,ij 

Ale 1 -* 



m 

— 7t teB In o 

8tt 2 s eBv 2 , 



+ 0(mi). (11) 



When the magnetic field is extremely strong, the wave- 
function renormalization behaves as = 8ir 2 m 2 /{eB) 
in the leading-order of eB, and the velocity is v 2 ^ — 
m 2 /{eB) ln(e _7 A 2 /m 2 ) accordingly. The contribution to 
the potential energy then becomes 



V v {m) 



dm' — - feB In 

8tt 2 s 



Ale 1 -"/ ln(e^A 2 /m' 2 ) 



eBm 2 

lD7T^ 



^Al 



(const.) , 



(12) 



where we keep only the dominant term for small m and 
drop negligible terms cx m 2 ln(lnm 2 ) and m 2 / lnm 2 . The 
potential contribution from the pion loops has a singular- 
ity at m = and leads to a divergingly positive curvature 
at small m, which favors chiral-symmetric phase in a way 
opposite to the quark potential in Eq. (5). Since V n (m) 
encompasses an opposite effect to the Magnetic Catalysis, 
we would call this the Magnetic Inhibition. Interestingly 
the logarithmic singularity associated with the Magnetic 
Inhibition is of the same strength as that with the Mag- 
netic Catalysis and these two effects should compete at 
sufficiently strong B. Our main purpose in this work 
is to propose a new physical mechanism, the Magnetic 
Inhibition, leading to a singularity ~ m 2 lnm 2 , and the 
determination of the singularity coefficient would require 
more works (and possibly depend on details of model as- 
sumptions). 



Here we would emphasize that our results are qualita- 
tively consistent with the latest lattice-QCD data indi- 
cating the decreasing behavior of chiral T c with increasing 
B. At T / the logarithmic singularity in Eq. (5) re- 
sponsible for the Magnetic Catalysis vanishes due to the 
absence of the Matsubara zero-mode for fermions and the 
Magnetic Catalysis is significantly weakened [20]. The 
Magnetic Inhibition is, on the other hand, enhanced by 
the temperature effects since the Matsubara zero-mode 
for bosons should be accompanied by a stronger infrared 
singularity. Therefore, the Magnetic Inhibition can soon 
overcome the Magnetic Catalysis at finite temperature. 
This could be a feasible explanation for the lattice-QCD 
data [8]. In fact, the lattice-QCD data implies that, for a 
fixed value of T, chiral symmetry is restored as a function 
of increasing eB. Quantitative analyses on this phase 
transition induced by the Magnetic Inhibition should de- 
serve more investigations including lattice-QCD simula- 
tions and chiral-model approaches. 

Finally, we would stress that the present work is the 
first attempt to exemplify the importance of the hadron 
structural change in a strong magnetic field. Our anal- 
ysis could be extended to investigate, for example, the 
possibility of the meson condensation induced by the B- 
effect [21]. 
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